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We show that the fractional statistics of quasiparticles in dimerized graphene, in recent proposals 
for charge and statistics fractionalization, can have two realizations depending whether elementary 
objects can be considered as point-like or extended objects. Therefore, there are two phases of 
proposed excitations and we give their topological descriptions with their respective statistics. We 
propose that a natural setting for fractionalization are certain systems with excitonic instabilities 
and demonstrate this by an example of graphene bilayer. 



I. INTRODUCTION 

Recent proposals^i^ for fractionalization in dimer- 
ized graphene and similar structures are based on the 
physics of charge fractionalization in polyacetylene~ But 
dimerization is hard to achieve in monolayer graphene. 
Mathematically fractionalization reduces to finding zero 
modes in a Dirac like equation for the dimerized pattern 
with vortex structure. A similar problem came from find- 
ing zero modes of vortex solutions of p-wave and other 
superconducting order parameters. From this comes the 
idea that in the systems needed for the fractionalization 
previously introduced, in which number of electrons is 
conserved quantity, excitonic instability with its BCS like 
hamiltonian and structure may induce order parameters 
and vortex solutions with only one zero mode necessary 
for fractionalization. Like quantum Hall bilayer at v = 1 
also graphene bilayer represents a natural setting for an 
excitonic instability as found in Ref. [||. And similarly 
but not obviously, as fractionalization in quantum Hall 
bilayer into merons with half electric charge and frac- 
tional quantum statistics, we expect similar quasiparti- 
cles - vortex solutions in graphene bilayer. 

In Section II we will review the recent proposals for 
fractionalization, where we will also clarify the situation 
and answer the question pertaining zero mode solutions 
and statistics that appeared in literature. Section III will 
be devoted to the bilayer graphene as a stage for charge 
fractionalization. 



II. FRACTIONALIZATION IN DIMERIZED 
GRAPHENE 

The effective (low-energy) Hamiltonian in the presence 
of the Kekule deformation of a graphene monolayer 
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u a ,v a and Ub,Vb denote the electronic (effective - Dirac) 
variables of the two triangular sublattices, a and b respec- 
tively, of the honeycomb graphene lattice. For the usual 
Kekule texture we get with A(r) = A G . With no tex- 
ture we have two cones in the spectrum, e±(p) = ±|p|, 
and with the Kekule texture mass gaps are opening : 
e ±(p) = ±-\/l-Pl 2 + |A | 2 in the single particle spectra. 



A. Charge fractionalization 

Let's assume a vortex structure in the complex param- 
eter A(r)i: 



A(r) = A(r) exp{-m0}, 



(2) 



where polar coordinates are used. We seek solutions for 
electronic states of the dimerized graphene in the pres- 
ence of this structure with zero energy. The equations 
that follow from Eq.([T]), in the case of sublattice a are 



d z u + iA(r)v = 
iA(f)u — d z v = 



(3) 



and with the exchange z «-> z they are also valid in the 
case of sublattice b. 

In the polar coordinates we have 

exp{— i9}(d r dg)u{r) + £exp{— in6}A(r)v(r) = 

r 

_ i 
ieKp{in6}A(r)u(r) — exp{i8}(d r H — dg)v(r) = 0(4) 

r 

In order to separate angular dependence we substitute 
u(f) — u exp{— wn#}tt(r) and v{f) = v exp{— U9}v(r) 
to have 

(d r Wr)+*A(r)— v(r) =0 

r u 



1 - V 

{dr- -)v(r)+iA(r)— u(r) = 
r v 
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if I = n — 1 — m. Further if we take iA(r) — — f(r) = 
|A(r)| that fixes the ratio ^f- the radial problem is re- 
duced to solving 



(Or 



(Or 

n - 1 



r 
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)u(r) + f(r)v(r) = 

)v{r) + f(r)u(r) = (6) 



There are two linearly independent solutions to the equa- 
tions. The behavior at large r is, apart from powers of 
r, exp{=F/^r}. Since the solution must be normalizable 
only one is acceptable. At the origin, the asymptots that 
follow, with /(r) ~ / rl n l, 



u(r) ^ Ul r m + u 2 r H+n - m 

n—l — m 



v(r) = Vl r^ +1+ 



v 2 r 



(7) 



To have single- valued, non-singular solutions at the origin 
we have to demand that m is integer and 



1 > m > 0. 



(8) 



The solutions that we get are similar but not the same as 
m Ref. 0. There superconducting couplings in the Dirac 
lagrangian induce different signs in the angular depen- 
dence u ~ exp{— im8},v ~ exp{i(rt — 1 — m)9} which 
then guarantees one angular momentum eigenstate per 
any odd value of the vorticity, v = —n. In our case only 
v = — 1 vorticity solution is an angular momentum eigen- 
value: m — 0. In this case the radial problem is simplified 
and the explicit solution is 



= C 





i exp{ia} 
1 




exp{- / f(r )dr } (9) 



where a is a constant defined by A(r) = 
| A(r)| exp{ia} — f(r) exp{ia} and C is a normal- 
ization constant. If we do not demand that the solutions 
have to be eigenstates of angular momentum the condi- 
tion, Eq.® ensures that we have n zero mode solutions 
in the case of the vortex with vorticity v = —n (n 
positive) 1 . Thus negative vorticity vortex states exist 
only on sublattice a, and very similar analysis shows that 
only positive vorticity vortex states exist on sublattice b 
in which case again there are as many zero modes as the 
value of vorticity isi. 

It can be argued^ that the charge bound to a vortex 
of vorticity v = 1 is — |. We have to study the change 
in the local density of states of the Dirac Hamiltonian 
Q] in the presence of mass twist in Eq. Q . Because of 
the sublattice symmetry^ to any negative eigenstate of 
the Dirac kernel ^ , _ e (r) corresponds to a positive energy 
state ty e (r), related to v I'_ e (r) by a unitary transforma- 
tion. Hence the local density of states 



v (j 



(10) 



is symmetric with respect to zero energy. Demanding 
the conservation of the total number of states after the 
inclusion of the mass twist, we get 



dr{2 / ^(r,e)de + |* (r)| 2 } = 0, (11) 



where ^(r) stands for the single zero mode. Its normal- 
ization to one leads to 



(12) 



df / 5v{r 1 e)de = 

J — oo 

so the net charge difference is — |. 

B. Statistics 



The introduced system can be in short described as 
Dirac electrons in the presence of a twisted mass^: 



C = y(ij u d v + Aexp{i7 5 0})* 



(13) 



where 7„, v = 0, 1, 2 are 4x4 Dirac matrices in the Weyl 
representation. The problem may be reformulated by 
dividing the vortex excitations into two groups, + and -, 
according to + and - value of vorticity corresponding to 
singularities in <p + and 0_ respecively, where cf> = cf> + + 
4>-, and introducing gauge fields^, 



a M = ^(<9 M + - ) 



(14) 



We denote by sgn{m^^ a quantity that at the singu- 
lar point of any vortex takes + or - depending whether 
the charge of the vortex is of sublattice a or b kind'. 
sgn(m$)Q^- represents vortex excitation electric charge 
current: 



r 



Ztt 



(15) 



sgn{m^) is necessary as we found out that vortex excita- 
tions with both positive (+1) and negative (-1) vorticity 
that live on sublattice b and sublattice a respectively, pos- 
sess the zero mode and unoccupied by electron represent 

— | absence of charge. 

sg^m^)^^- represents current of axial charge (valley 
index) 2 associated with vortex excitations (normalized by 
vortex axial charge) - j£. The sign of axial charge in the 
density-current j§ comes from sgn{m^,). This comes from 
the correspondence - see Appendix, in the sign of the sub- 
lattice density difference 'J 73 ^ (our states are eigenstates 
01 7o73-) and the sign of the expectation value of ^7075^1' 

- axial charge density. The expectation value of the axial 
charge is ±i like of ordinary charge. 

Then a simple statement follows for the topological 
part of the effective action in the dual representation of 
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the theory i.e. in terms of vortices instead of Dirac par- 
ticles (electrons). (In the dual picture elementary (2ir) 
fluxes of gauge fields are particles and the gauge fields 
represent particle background.) The form of the topo- 
logical part of the action is 



sgn(m 3 ) 1 - 1 

-A_ a(d x b) - -af - - 



(16) 



We just encoded in Eq.(fTB]) the expressions of our cur- 
rents found previously. There is overall \ factor because 
of the value of the charge of the vortices that must cou- 
ple as — ^A^j^ to the external field if we introduce it 
(a -> a + A). 

If we introduce gauge fields 

R = a + b,L = a-b (17) 
we arrive at the following form of the Lagrangian, 



sgn(m 3 ) 



,3~ 



if 



8tt 



where 



and 



(RdR-LdL)-R( J Y )-L("- 2 



-■+ 3 + 3s 
J 2 
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(18) 



(19) 
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the currents of the good quantum number - vorticity as 
opposed to charge j and axial j'5 currents. Considering 
the Aharonov-Bohm phases for encircling quasiparticles 
around each other we easily and clearly get that quasipar- 
ticles, j + and j~ , have semionic statistics among them- 
selves and their mutual statistics is trivial. 

Therefore through a straightforward analysis of quasi- 
particles as point particles i.e. singularities of phase - 
vortex solutions we come to the conclusion that they 
obey semionic statistics an d the ir theory is the doubled 
Chern-Simons i.e. t/(l)2 x U(l) 2 - see Ref. [HI as found 
in Ref. 0. On the other hand in Ref. the quasiparticles 
were viewed as extended objects - meron configurations 
of field 71, n 2 = 1, and the conclusion was that they pos- 
sess quarton statistics. If we apply a simplification that 
vector n is always in the x — y plane except at the center 
of excitation we will find following the arguments of Ref. 
that the excitation possess semionic statistics. Consid- 
ering this we may ask ourselves what modifications of our 
approach are necessary to account for extended vortices. 

The electric charge current j M and axial charge current 
jg can be introduced in the topological part of the action 
by simply taking j M = ij M and j'5 = hj^ because the 
point charges carry half of the unit of electric and axial 
charge: 



sgn(m 3 ] 

sgn(m 3 ] 
8tt 



a(8 x b) - a^f - 
(RdR - LdL) - Rj + - 



Lj~. (21) 



Now we can see what principle can guide us to modify the 
theory. j + and should, in principle, correspond to the 
charged fermions, electrons that may appear even in low- 
energy theory (they should certainly exist in a complete 
theory). To have that we will add appropriate Chern- 
Simons term as additional dynamics that is allowed : 

S9nM -a(d x b) + S J^a(d x b) - a,f - b^. (22) 



2tt 



2tt 



Now fractionalized excitations, 
have quarton statistics i.e. 



2j» and 



2/ 
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S J^ ( R d R-LdL)- 1 -Rr- 1 -L]-. (23) 

This theory alone (as a topological one - so-called BF 
Chern-Simons field theory and without sgn(m 3 ) which 
can be absorbed by simple redefinitions) was investigated 
in Ref. |ll| as the theory of 2d superconductors with vor- 
tices and quasiparticles as excitations, and was proposed 
as the description of the topological part of a phase for 
the QH bilayeri 12 These considerations also imply that 
quartons can be only found confined in pairs. 

The beginning Lagrangian in (|13|) can be restated by 
a gauge transformation: &± — * e z ^ >± ^±, where &± arc 
chiral components of Dirac field, ^± = ^(1 ± 75)^ as 2 



£ = tf(i JB- A -75 ^ + A)*. 



(24) 



After integrating out Dirac fermions, the total C can be 
expressed also as 2 



C 



12A 
1- 



(d x a) 2 
1~, 



2tt 



sgn(m 3 ) 
2^ 



a{d x b) 



^f a »-\3t h n- 



(25) 



When A or screening charge (in the Maxwell term) is 
large we may expect that the point-like description (via 
semions) of vortices is appropriate, but when A is small 
we are in a superfluid phase where presumably quartons, 
but certainly some extended objects, are confined and 
appear in pairs. This latter physics remind us of the 
quantum Hall bilayer physics where merons in the super- 
fluid phase for the bilayer are extended objects, appear 
in pairs, and have quarton statistics^. And indeed the 
action (after integrating out fermions in the presence of 
an additional field - a staggered chemical potential i.e. 
a third component of the ft vector) that was found in 
Ref. 0, is similar to the effective action for the quan- 
tum Hall bilayer (as in Ref. [l3T ) where also the fermionic 
current is equal to the 0(3) topological current^. The 
theory in Eq. (|22[) maybe a crude oversimplification, but 
it tells us what is the feature of any phase which includes 
both, fractionalized excitations (merons) and fermions - 
the excitations must be bound in pairs (compare Ref. 0). 
In other words when merons are deconfined, they can 
be viewed as point-like objects with semionic statistics. 
That was also a result of the numerical study in Ref. [U 
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of small (quantum dot) systems with spin in which de- 
confincmcnt of merons was proposed and mapped to a 
spinon (semion) gas of Haldane-Shastry chain. 

The conclusion is that the difference between Ref. 
and Ref. in assigning the exchange statistics comes as 
a difference in how one considers n: (1) as a continuous 
vector field as in the quantum Hall bilayer where quar- 
ton statistics (meron - extended description) will survive 
even with no bias that is, in Ref. 0, with no uniform 
fi s - staggered chemical potential but with adjustments 
of /j, s at places of excitations, or (2) as a field that has 
to take values only in the x — y plane as in Ref. [2j (and 
only uniform, constant value of /i s is allowed) with sin- 
gular behavior on vortices and hence scmionic statistics. 
So the question is whether we are promoting fj, s into a 
dynamical variable. Both possibilities seem allowed but 
lead to different phases in general, in which fractional 
objects have same charge but different statistics, and are 
confined and deconfined respectively. 

The derivation of the semionic statistics (Eqs. [14] - 
UnUSO])) that we gave is very general and still valid even 
with inclusion of a time reversal breaking term - n 5" 75^ 
that was included and discussed in Ref. 0. The only 
assumption is the value of the charge of the excitation 
for that case, which we can safely take to be h in ac- 
cordance with the n formalism in Ref. (for )a s = 0). 
A related question or comment may be that the gauge 
transformation of Ref. of C in Eq. (TIB"]) will not lead to 
a simple transformed form with only fields: ^ , a M , and 
(see Eq. ([24|) ) in this case. That is true but it is very un- 
likely that (small) 77 perturbation will lead to an effective 
Lagrangian of ^>,a^, and fields with a change of the 
coefficient of the minimal coupling of field a M from one to 
two that is needed for quarton statistics if we follow the 
same steps as in the derivation of the semionic statistics. 
Therefore, the time reversal breaking term alone can not 
lead to a statistical transmutation of semions into quar- 
tons, although it seems its presence is the only way to 
recover and demonstrate (Ref. 0) quarton statistics in 
the ft formalism ( when (A > |?7|)). The inclusion of 
the uniform staggered chemical potential will change the 
overall factor ^ (the charge of the excitation) that multi- 
plies E ci. (TIT)]) , but still the statistical angle will be given 
by - = Q s where Q s is the charge of the excitation ( and 
notby I = Ql). 



III. FRACTIONALIZATION IN BILAYER 
GRAPHENE 

The question is what are physical systems that may 
support fractionalization - obviously dimerized graphene 
is a hypothetical system. The bosonic degrees of freedom 
that we need can come as a result of electronic correla- 
tions, most notably excitonic and superconducting. From 
these, only excitonic conserve charge and may produce in 
their defects charged zero modes (as opposed to neutral 
zero modes in superconductors) that we need to have 



fractionalization. 

In order to facilitate the discussion of the excitonic in- 
stability and its zero modes in graphene bilayer, we will 
first discuss zero modes in the case of p-wave supercon- 
ducting and excitonic system. 

(1) The effective BCS Hamiltonian for the quasiparti- 
cles is 



k 



H = y>c+ Cfc + x(A*c_ feCfe + A k c+c+_ k ). (26) 



Introducing the Bogoliubov transformation: 



ctfc = u k c k - v k c_ k 



(27) 



that should diagonalize the Hamiltonian into H = 
Y] E k a\a k + const, implies Bogoliubov-de Gennes equa- 
tions 



E k u k = £ k u k - A* k v k , 
E k v k = -£, k v k - A k u k . 



(28) 
(29) 



If A k — A(k x — iky) and £ k « — /i(/i > 0), in the long- 
distance approximation the equations for zero mode(s) 
become: 



[iu — A{—i)c\v = 0, 
uv — (—i)Ad z u = 0. 



(30) 
(31) 



We may ask for zero modes that exist in vortex so- 
lutions for which we demand v(9 + 2n) = —v(9) and 
u(8 + 2tt) = —u{ff) in polar coordinates^. By solving 
(150)1 with A = const we are neglecting the short (small 
radial) distance details of the solution. We seek the so- 
lution in the following form 



u(r) v(r) 
and the equations that we get are 



^^-+iA^ T e l6 d r v{r) = 0, 
+ iA^e' l9 d r v(r) = 0. 



(32) 



(33) 
(34) 



To separate angular and radial dependence we must 
have k = I = i and the equations are reduced to 



fiu(r) + iAd r v(r) = 0, 
/iv(r) + iAd r u(r) = 0. 



(35) 
(36) 



If u(r) — uof(r) and v(r) = vof(r), the equations 
reduce to a single one, 
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Hf(r) + Ad r f(r) = 0, 



(37) 



if Ha — _ j_ Therefore our solutions can be cast in the 



following form: 



le 4 — =-, v 



(38) 



3 r 



where /(r) is of the simple radial dependence ■~ 
for /i = const. 

The usual approac h 16 ' 17 is to model order parame- 
ter with the vortex singularity i.e. to take, instead of 
Afe, e 4 2 Afce l 2 in our case. This symmetrized with re- 
spect to phase expression is used to ensure the anti- 
symmetry of the order parameter i.e. that the term 
J dr4 r+ (r)AfeW''(f) in the Bogoliubov Hamiltonian is 
well-defined and consistent with the anticommutativity 
of fermi operators when phase is coordinate-dependent. 

From the approach we used in getting the zero mode 
we can turn to the usual approach by the simple phase 
transformation u — > e~'5ii, v — > e l ^v, so that at the end 
our zero mode solution has the components: 



. ilt /(r) iJt /(r) 
u = ie l ? J -A-L, v = e l ~ ±AJ. 



(39) 



Now the quasiparticle operator for the zero energy 
state can be written as 



7o~ = / d 2 r(u(r)c + (r) + v(r)c(r)) 



(40) 



and immediately we can conclude that 7q = 70 for our 
solution i.e. it represents neutral Majorana mode. 

We should notice that with respect to the Dirac prob- 
lem in the dimerized graphene here, in the latter ap- 
proach, momentum operators are together, in the same 
term, with order parameter and phase singularity in the 
Hamiltonian. 

On the other hand, in the case of the excitonic problem 
which may be related and become a physical realization 
of dimerized graphene lattice^, we have a different basic 
Hamiltonian: 



k 

Introducing 



/■ 'Ik Ik) 



£(A fe /?+ 7fe + Afr+flO 

k 

(41) 



e k u k = E k u k + v k A* k , 
e k v k = E k v k -u k A k . 



(44) 
(45) 



For the zero modes (e* = 0) simple redefinition v k — > 
—v k transforms the equations into the same as for the 
superconducting problem. 

We will assume A k = A(k x — ik y ) and E k = e = 
const > 0. Then for the zero modes we have 



(46) 
(47) 



eu + A(-i)d?v = 0, 
ev — A(—i)d z u = 0. 



The equations are the same as in Eq. ([30j) and Eq. ([3Tj) . 
Again, we may ask for the zero modes that exist in vortex 
solutions for which we demand v(9 + 2ir) = —v(0) and 
u{6 + 2tt) = —u{&) in polar coordinates and the answer 
would be the same. But if we want to stay in the language 
of the order parameter, we may model it, in the presence 
of a vortex, as A k c l8 , with no symmetrization as was 
necessary in the superconducting problem. In this case 
the solutions become 



f(r) 



(48) 



which are not satisfactory for the electronic wave func- 
tions because they are not single-valued while the singu- 
larity is borne by the order parameter. If we choose in 
the order parameter e~ 1 ^ A k e 1 ^ , our solutions are repre- 
sented by 



u = ie^e ie ^i,v 



(49) 



and this, although not an angular momentum eigen- 
state, is the solution that can describe a charged mode 
because u* / « as a crucial difference with respect to 
the superconducting case. The excitonic Bogoliubov 
transformation mixes the same kind of charged opera- 
tors, leading to charged zero modes as in the dimerized 
graphene problem. Although the precise form of the 
equations for the zero modes is not the same in the two 
cases (excitonic p-wave system and dimerized graphene), 
they are very similar in the long-distance regime and 
should lead to the same conclusion about the charge frac- 
tionalization. 

(2) For the case of bilayer graphene excitonic insta- 
bility described in Ref. d, we have in the limit of small 
interlayer hopping 



B k = u k f3 k - v k j k , 
C k = v k p k + u k j k 



(42) 
(43) 



that should diagonalize the Hamiltonian into H 
E k tk(B+B k -C+C k ), we get 



2^i 



(50) 



for the energy in Eqs.(|44p. Here t± is the interlayer 
hopping parameter, ±V is the bias that causes excess 
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electrons and holes in the lower and upper graphene lay- 
ers, respectively, and is the bare kinetic energy for 
which Ek ~ k. The gap function is of the following form, 
near two nodal points^, 

A fc = i(k x t ik y )\k\A, (51) 

where A is a positive constant. We will write as 
Ek = V —5k 2 with also 5 being a positive constant. Then 
the equations (|44|) with = become 

To consider the vortex solution we have to set 
e"'5Aie'5 instead of Afc, which if we want to absorb 
the phase into v and u leads to slightly different equa- 
tions due to the presence of the k 2 terms that we also 
must keep in this long-distance analysis. Because 



We will assume the asymptotic behavior of u(r) and 
v(r) as u{r) ~ e~ Ar and v(r) ~ e~ Ar , and justify it in 
the end. Then, for example, in the case of v(r) we have 

I Q i . I d d 1 z 1 / 2 

vtzw** drV{r) = - Xvo y Tzdk { **n 

oo 1 

X E ^(-A)"z"/ 2 z" /2 } 

n=0 

1 / d I 8 

= - Xv °7^V z dzVaz- (59) 

1 °° 1 

V^i ^ Q n! 
In + l ~ In — 1 i+i_fc_i^i ,„„,. 

x \/ — - k \l— z 2 * 2 • (6i) 

Therefore the Eq. ([55|) becomes 



e 2 — e 

9z <9z 



= e 2 — e 2 e 2 — e 2 

= r— -— )(— —) 



the equations for the vortex solution become 



., d l.,9 1 ' 



(54) 



(55) 



(56) 



V + S(- 



2z'^dz + 2z> 



u(r) 



+Az(-AH-^E^(-A) 



ZZ — 

Tl — 



71 + 1 



-2 2 K Z 2 + 



(62) 



To factor the angular dependence we choose k = 1 = j 
and in the long distance approximation, for which we take 
also \J n(n — 1) ~ n in the series expansion, we have 



[V + 5d 2 ] u + iA(-\) 2 v(r) =0, (63) 
[^ + 5<9 2 ] u - iA(-A) 2 u(r) =0. (64) 

This means 



with the requirement u(8 + 2tt) = —u{6) and v(# + 
2ir) = —v{6). If we seek solutions in the following form 
(see Eq.jSSJl) 



u(r) v(r) 

1-,V = 7- 



(57) 



we face the problem of properly performing the oper- 
ations with the square root operator, \J ■§£■§=■ Its action 
we will define on the space of monomials in z and z, for 
which we have 



dz V dz 



(58) 



mo [V + <5A 2 ] + iA\ 2 v a = 0, 
v [V + 6X 2 ] - iAX 2 u = 



A? 



V \2 V 

-(8+A) ° r A 2 — ~A=g- 



(65) 
(66) 
(67) 

Both, Ai = ±i^Jg^ 5 , 

represent delocalized zero modes. In the case of A2 that 
would be if S > A, otherwise we have only one physical 

localized state with e~ A2r ,A2 = \J , decay function. 

Comparing with solutions in Refill for A for fixed pa- 
rameters for the graphene bilayer, we find that for large 
enough V , bias parameter, we can have the bound solu- 
tion. Then the three zero modes carry fractional charge, 
I, that we get applying the same arguments that were 
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given m RefQ 

for the dimerized graphene. Essential in 
these arguments is that Dirac kernel has the correspon- 
dence of positive and negative energy solutions. In this 
case we also have the Dirac structure of the problem, 
with the nodal points, and the same correspondence. 

In this Hartree-Focb^ treatment of the excitonic prob- 
lem we get three zero modes (for large enough V) and 
the question comes whether they will stay with fur- 
ther inclusion of interactions. In the case of three zero 
modes, the charge that they carry is mostly smeared out 
through the system. If the interactions were able to split 
the two delocalized zero modes, we could have localized 
charge | vortex excitation similarly to the case of dimer- 
ized graphene. 1 The conclusions about the topological 
terms and fractional statistics that we reached in the 
case of dimerized graphene will still hold. The interaction 
term that we have in mind would be ~ V pk o p-k a where 
k = | Ai| and an excitonic coupling between the mode k a 
and — k Q would produce necessary splitting. Otherwise 
(with no splitting) the situation is less clear but if we 
assume that the localized mode describes a missing of | 
charge the delocalized modes will describe additional de- 
grees of freedom (that may be occupied or unoccupied) 
that may further decrease the value of the statistical an- 
gle of excitations or equivalently increase the coupling of 
the doubled Abelian Chern-Simons term. 

The recent work, Ref. [H, that appeared while we 
were finishing the writing, concerns fractionalization in 
excitonic bilayer graphene that is not naturally (Bernal) 
stacked but consists of two parallel layers at some larger 
distance that leads to even number of zero modes due to 
the valley degeneracy and, therefore, no charge fraction- 
alization. Still we can not rule out, on the basis of the 
long distance analysis of vortex solutions in the excitonic 
condensate that we presented, the same doubling in our 
case. This important question can be resolved only by a 
detailed, numerical analysis of the bilayer graphene. 



where we simplified the decay function by taking | A(r) | = 
|A| = const, and again the C is the normalization con- 
stant. Notice the absence of two different operators (one 
for particle, the other for hole) that we would have for 
a non-zero energy level. For the sake of the argument, 
going in reverse, we can fix the normalization constant 
by demanding that the charge associated with the zero 
mode is i- e - 



<fr < *+(f)W (r) >= 



(A.2) 



That would imply 



dV<$ +(rV* (r>=-, 



(A.3) 



for the expectation value of the axial charge, where we 
took 



75 



10 
10 
0-10 
-1 



in the Weyl representation. On the other hand our state 
is an cigenstate of the sublattice charge difference oper- 
ator: 



R = "3 = 7073 



IV. ACKNOWLEDGMENT 

The author thanks Chang- Yu Hou for correspondence. 
The work was supported by Grant No. 141035 of the 
Ministry of Science of the Republic of Serbia. 



" 





1 


0" 




" 





-1 


o - 




" 1 




















1 













1 







-1 








1 













1 



















-1 





. 


1 





. 




. 


-1 





. 




. 








1 



APPENDIX 

In the second quantized formalism the zero mode so- 
lution in Eq. ^ , of definite vorticity v = — 1, contributes 
to the expansion of the Dirac field as a term equal to 












"0" 








i exp{ict} 



co + 
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4 > 











.0. 





(A.l) 



so that 

J d 2 r< *+(r)a 3 * (r) >= 1 (A.4) 

of the same sign as the expectation value of the axial 
charge. Both signs would reverse if the solution were 
with vorticity v = 1, on sublattice b, although the sign 
of the electric charge would remain the same. 
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